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ABSTRACT 

We have used a new observing system on the Parkes radio telescope to carry out a series of pulsar observations 
of the globular cluster 47 Tucanae at 20-cm wavelength. We detected all 1 1 previously known pulsars, and have 
discovered nine others, all of which are millisecond pulsars in binary systems. We have searched the data for 
relatively short orbital period systems, and found one pulsar with an orbital period of 96 min, the shortest of any 
known radio pulsar. 

The increased rate of detections with the new system resulted in improved estimates of the flux density of the 
previously known pulsars, determination of the orbital parameters of one of them, and a coherent timing solution 
for another one. Five of the pulsars now known in 47 Tucanae have orbital periods of a few hours and implied 
companion masses of only ~ 0.03 M0. Two of these are eclipsed at some orbital phases, while three are seen at 
all phases at 20 cm but not always at lower frequencies. Four and possibly six of the other binary systems have 
longer orbital periods and companion masses ^ 0.2 M0, with at least two of them having relatively large orbital 
eccentricities. All 20 pulsars have rotation periods in the range 2-8 ms. 

Subject headings: binaries: eclipsing — binaries: general — globular clusters: individual (47 Tucanae) — 
pulsars: general 



1. INTRODUCTION 

Soon after the discovery of the first millisecond pulsar 
(Backer etal. 1982), Alparet al. (1982) proposed that low-mass 



X-ray binaries (LMXBs) are the progenitors of millisecond pul- 
sars. Given the high incidence of LMXBs in globular clusters, 
searching for millisecond pulsars in clusters became a natural 
goal. After much effort, the first such object was discovered, a n 
isolated pulsar with a period of 3 ms in M28 (Lyne et al. 1987). 
Following searches of many clusters in the late 1980s and early 
1990s, over 30 pulsars are now known in 13 clusters. In con- 
trast to most pulsars in the Galactic disk, the majority of cluster 
pulsars have short rotation periods and are ofte n members of 
binary systems (see Kulkarni & Anderson 1996 for a review). 



Searches for millisecond pulsars in globular clusters en- 
counter two main difficulties. Firstly, since globular clusters 
are on average more distant than the ^ 2 kpc within which most 
large-area surveys have been sensitive t o millisecond pulsars 
(see, e.g., Camilo, Nice, & Taylor 1996), much longer obser- 



vation times are typically required than for Galactic searches. 
This in turn increases the complexity of the data-reduction task. 
Such observation times are feasible only because of the small 
number of search targets. Secondly, for binary pulsars there is 
the added difficulty that the apparent pulse period may change 
significantly as the pulsar experiences orbital accelerations dur- 
ing the longer integration. If this issue is not addressed at the 
data reduction s tage — and with significant exceptions ( Ander- 
son et al. 1990) it rarely is, because of the massive compu- 
tational requirements — many relatively bright binary pulsars 



may remain undetected. 

Despite these difficulties, the rewards of a thorough search 
can be significant. Because of the high probability of close stel- 
lar interactions, unusual binary systems not otherwise expected 
in the Galactic disk may form in clusters. In addition, for clus- 
ters where several pulsars are found and studied with the detail 
allowed by high-precision timing measurements, much can be 
learned about the clusters themselves. Applications to date in- 
clude a determination of the core mass distribution from mea- 
surements of pulsar accelerations and positions, and limits on 
intra-cluster me dium propertie s , by measuring t he dispersion of 
various pulsars (Phinney 1992; Anderson 1992). 



The basic observed and derived properties of 47 Tucanae 
(NGC 104) used throughout this paper are summarized in Ta- 
ble [l| The fact that it is a relatively nearby cluster with a high 
stellar density makes it an attractive target for pulsar searches. 
47 Tuc together with M15 have been by far the most prolific 
cluster targets for pulsar searches to d ate, with surveys reveal 
ing 11 and eight pulsa r s respectively ( [ Manchester et al. 1991 
Robinson et al. 1995; ^nderson 1992 ). The observed popu 



lation of pulsars in these two clusters is very different: while 
those in 47 Tuc are all millisecond pulsars, the pulsars in Ml 5 
are more varied in their properties, including slowly-spinning 
puls ars and a rare double neutron-star binary (P rince et al. 
1991). It is also worth noting that although all pulsars in M15 
have had positions and period derivatives determined using pul- 
sar timing techniques, this has been possible to date for only 
two of the pulsars in 47 Tuc. This is due mainly to the different 
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Table 1 

Observed and derived properties of 47 Tucanae. 



Reference 



Center R.A., a„ (J2000) 
Center Decl., 547 t„c (J2000) 

Distance 

Central mass density, p(0) . . 

Angular tidal radius^ 

Angular core radius"^ 

Escape velocity 

Velocity dispersion 



00''24"'05?9 
-72°04'51"1 
4.6 kpc 

4x lO^M0pc"- 
40' (54 pc) 
12;'2(0.27pc) 
56.8 km s"' 
13.2kms"' 



h„t,„tt,aV„rtn ,t „1 1QQT 






ftal. 199/ 


|webbink 1985 













Dp Marclii pt ^ 


1. 199( 








Webblnk 1985 



*The value quoted is a revision of Webbink's (1985) catalog entry (see for details). 

^ Also known as the "limiting radius", this is the radius at which the King model surface brightness-radius profile drops to zero ( King 1962 1. 

"^The core radius is a scale factor used in the King model, and denotes the radius at which the number density of stars in the cluster per unit area is half the central value. 



telescopes used and different interstellar scintillation character- 
istics: at Arecibo, Puerto Rico, the majority of the distant pul- 
sars in M15, whose apparent fluxes vary little, are detectable on 
most days, while at Parkes, Australia, most pulsars in 47 Tuc 
are visible only when the interstellar weather is propitious and 
the pulsar signals are amplified in the observing radio band by 
scintillation. Unfortunately this has meant that many of the sci- 
entific objectives suggested by the wealth of pulsars known in 
47 Tuc remain unfulfilled. 

In mid- 1997 a new, sensitive system for observing pulsars at 
a wavelength of 20 cm became available at the Parkes telescope. 
Puls ars are in general steep-spectrum sources (L orimer et al. 
1995), and all pulsars in 47 Tuc were discovered at a wave- 
length of 5 or 70cm (Manchester et al. 1990, 1991; R obinson 
et al. 1995). However, the new system provided such perfor- 
mance improvements that observations to study the known pul- 
sars and to find previously undetected pulsars seemed desirable. 
These observations were so successful in the rate of detection 
of previously known pulsars and with discovery of new ones, 
that we have moved to observing the cluster almost exclusively 
with the new system. In this paper we report on the discovery 
of nine binary millisecond pulsars in 47 Tuc and give improved 
parameters for some of the previously known pulsars. 

In §^ we describe the observing system used, list the obser- 
vations made, and outline the data reduction carried out. In §|^ 
we estimate the sensitivity of the survey for isolated and long 
orbital period pulsars, and compare it to that of past surveys 
of 47 Tuc. In §^ we summarize the results, including param- 
eters for the newly discovered pulsars and updated parameters 
for some previously known systems. Finally, in §^ we discuss 
some limitations, and implications, of our survey and its results. 

2. OBSERVATIONS AND DATA ANALYSIS 

Since 1997 October we have been using the Parkes tele- 
scope together with the pulsar multibeam data-acquisition sys- 
tem to observe the globular cluster 47 Tuc at a central radio 
frequency of 1374 MHz. The system consists of a cooled, dual 
linear pol arization, 13 -beam receiver package ( Staveley-Smith 
et al. 1996), and new filter banks, single-bit digitizer, and data- 
acquisition computer 

Each of the 13 beams has a half-power diameter of about 14' 
on the sky. Since a single beam encompasses the great majority 



of the cluster stars (cf. Table [i|), we record data from the cen- 
tral beam only (but see discussion in § 5 . 1 .2 ), while pointing the 
telescope at 47 Tuc C, the first pulsar discovered in this clus- 
ter. Total-power signals from each of 96 contiguous frequency 
channels, each 3 MHz wide, are one-bit sampled every 125 /is 
and, together with relevant telescope information, are written 
to magnetic tape for off-line processing. Data are usually col- 
lected on d ays when a pulsar survey of the Galactic plane is also 
underway ( Lyne et al. 2000[ ), and are recorded from 47 Tuc for 
up to five hours daily. Table || contains a log of all observations 
reported in this paper 

The largest continuous data sets recorded are 4.66 h long, 
corresponding to 2^^ time samples. With the computing re- 
sources presently at our disposal, it is impracticable to reduce 
such data sets in single, coherent, Fourier analyses. Further- 
more, since about one third of the pulsars already known in 
47 Tuc are in binary systems, s ome with orbital periods of un- 
der 3h ( Robinson et al. 1995[ ), we would select against such 
systems if we analyzed each large data set as a whole. As a 
compromise between computational feasibility and sensitivity 
(both to weak radio sources and relatively short-period bina- 
ries), we restrict the data analysis reported in this paper to con- 
tiguous blocks of 2^-' time samples (17.5 min). 

The data were de-dispersed to a reference dispersion mea- 
sure (DM) of 24.6cm-^pc, the DM of 47 Tuc C. This was 
done by delaying samples from the higher-frequency channels 
to account for the faster propagation through the ionized inter- 
stellar medium, before summing data from all channels into a 
one-dimensional time series. For the present analysis, where 
we are particularly interested in looking for binary pulsars with 
short orbital periods, this choice was made largely to keep the 
total computational effort manageable. The pulsars presently 
known in 47 Tuc exhibit a range of DMs between 24.1 and 
24.7 cm"^ pc (Table |[). In the current analysis, any pulsar with 
a true DM that is within one unit of 24.6cm"^pc will have its 
pulse smeared by at most I ms as a result of this single-DM 
approximation. 

A s noted by a number of authors (M iddleditch & Priedhorsky 
1986; pohnston & Kulkarni 199 1| ), standard pulsar search anal- 
yses, which look for significant harmonics in the power spec- 
trum of a de-dispersed time series, suffer reductions in sensi- 
tivity to pulsars in short-period binary systems. The effect of 
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Table 2 

Log of observations and detections. 





Dale 


Number t)f 


c 


D 


E 


F 


G 


H 


I 


J 


L 


M 


N 


o 


Q u 


(MJD) 


iLiiL-graiioLis 




























50683 


16 


0/4.6 


2/7.0 


14/19.8 


0/2.3 


0/2.5 


10/17.0 


0/2.4 


6/8.8 




0/2.5 




0/2.9 




50686 


16 


4/7.6 


0/3.2 


0/4.4 


16/33.5 


0/3.4 






16/12.2 


0/3.1 






4/10.1 




50689 


16 


5/9.8 




0/3.1 


0/3.4 


0/2.9 




3/6.0 


11/12.9 


0/2.8 




16/18.2 


0/2,7 




50690 


16 


16/21.9 


0/2.8 


0/3.2 


0/6.1 


0/3.3 






1/6.6 








3/8.6 


0/2,0 


50739 


5 


5/20.9 


3/10.2 










0/7.4 


5/31.5 












50740 


8 


0/4.4 


1/7.2 


0/7.2 




0/3.7 




2/8.8 


8/13.4 




0/3.8 




0/5.3 


7/13,5 0/4,2 


50741 


8 




8/11.9 


8/19.7 








2/9.8 


8/46.5 








0/4.1 


0/4.5 0/3,0 


50742 


8 


8/15.7 


0/4.1 




2/6.5 


0/3.1 


8/16.6 


0/6.1 


6/12.6 


0/4.0 








0/4,0 


50743 


3 


0/7.9 






3/21.9 


3/12.8 






3/17.8 












50744 


8 


8/15.3 




0/3.2 


0/3.6 


8/14.7 


0/4.2 


0/2.9 


4/11.1 








0/3,2 




50745 


8 


0/2.8 


8/22.6 




0/5.4 




5/11.3 


0/3.7 


8/26.7 








0/3,5 




50746 


8 


0/5.4 


8/17.7 




0/3.1 




0/6.5 


0/3.6 


8/32.2 










0/6,0 


50748 


4 


4/10.5 


0/3.7 






0/3.9 










4/14.5 




0/4,3 




50980 


8 


8/22.5 


8/33.4 


8/20.7 


1/4.9 








0/3.8 






0/3.4 




0/2,7 1/8,4 


50981 


16 


16/63.2 


5/9.3 




15/18.8 


1/4.2 






7/7.9 




1/6.3 




0/3,0 




50982 


16 


2/5.8 


16/22.7 


9/13.3 


1/5.7 




0/2.5 




3/5.4 




1/8.5 




0/4,3 




50983 


4 


0/5.0 




3/12.6 


0/6.3 
















0/7,4 


0/5,6 


50992 


8 


0/5.4 


0/6.7 




6/12.9 








0/4.7 










0/3,3 


50993 


8 


8/40.2 


0/7.2 


0/5.4 


0/2.9 






1/9.1 


8/22.8 


7/11.8 








0/3,0 


50998 


16 


0/2.1 


3/6.6 


1/5.4 


2/5.4 








0/2.7 


0/2.7 




14/11.8 






50999 


5 


0/6.4 


5/14,6 




4/10.8 






0/4.1 


0/6.8 








0/4,8 


0/3,7 


51000 


16 


0/6.7 


0/2.1 










0/2.7 


2/6.7 


0/3.1 




0/2.3 


0/2,6 


0/2.2 


51001 


10 


5/10.2 


0/3.2 








0/6.0 




0/3.1 










1/2,9 4/10,9 


51002 


16 




2/5.5 


1/7.4 








0/2.8 


16/146.3 




0/2.2 




3/6.4 


3/7,7 


51003 


16 




4/8.6 


4/8.8 


0/6.1 


0/2.8 


0/3.2 


0/2.6 


16/36.1 








0/4.9 


0/5,3 


51004 


16 


0/5.8 




0/5.7 


12/12.6 








13/16.2 










10/10,2 0/3,0 


51005 


12 


2/7.7 


0/3.0 




0/2.5 








9/12.8 


0/2.7 






2/4.4 




51007 


10 


0/4.3 


0/5.2 












10/23.5 










0/2,9 


51012 


8 


5/11.4 




0/7.1 


4/10.9 


1/6.2 






8/186.1 








3/7.2 


0/4,7 


51026 


16 


16/25.1 


0/4.4 


0/4.0 


0/3.0 


0/3.6 






4/8.3 




0/2.2 




0/3.2 


0/2,4 


51027 


7 


6/15.7 


7/36.0 


3/9.9 








0/7.2 


7/14.3 




0/3.4 






0/3,1 


51029 


5 


0/4.4 




1/9.5 


0/5.4 






0/4.1 


5/26.2 












51031 


6 


0/3.1 




0/4.7 


4/10.2 






0/6.7 


3/8.0 








0/5.1 


0/7,1 


51033 


6 




3/12.8 






0/5.5 






0/3.4 










6/20,7 


51036 


7 


0/3.5 


7/15.2 




4/9.9 








0/4.5 








0/5.8 


2/9,6 0/3,9 


51037 


16 




16/34.1 








0/2.6 


0/4.6 


16/29.7 




0/3.8 






0/2,2 


51038 


6 


6/18.3 


5/10.2 




1/8.1 






0/3.8 










0/6,2 


10/10,7 


51039 


3 


3/17.2 


3/19.9 


1/8.6 






0/8.4 


1/10.3 














51040 


4 


4/18.0 


4/16.8 






2/12.1 




0/8.3 


4/42.0 








0/5,0 


0/8,2 


51085 


16 


5/8.9 


10/9.6 


0/2.8 


0/3.0 






0/2.1 


16/16.5 




0/3.2 


0/3.2 


0/4,4 


0/3,4 
































51088 


19 


0/2.2 


0/2.9 




0/2.8 


0/1.9 


1/5.7 


0/1.9 


8/13.0 




0/2.2 




0/2,7 




51089 


4 


4/13.6 






0/4.8 






3/10.3 


0/5.1 




0/4.6 








51090 


4 


4/23.2 


4/20.5 




0/5.7 








4/16.7 










0/7,9 


51091 


4 


4/49.3 


4/21.3 








4/21.2 




1/9.3 












51092 


8 


8/29.4 


7/15.5 


1/7.7 


8/29.3 






5/13.4 


0/4.8 




0/6.4 






0/2,8 0/4,5 


51093 


8 


0/6.8 


0/3.5 




0/3.2 






0/4.6 










3/6.7 




51095 


4 


4/18.0 




0/6.8 








3/11.8 


1/10.7 












51096 


5 


4/30.5 


5/16.4 




0/4.5 








0/4.8 












51097 


12 


3/9.8 


4/12.2 




7/18.2 




1/5.2 




12/88.3 




8/11.2 




0/3.7 




51098 


16 


15/17.8 


0/2.4 


0/3.7 










0/2.9 


2/6.6 


2/6.2 








51099 


19 


0/2.7 


4/6.9 


18/28.6 


1/3.8 








0/3.6 








0/2.6 


10/8,2 


51100 


16 


0/3.5 




0/2.4 




0/3.5 


1/6.0 




16/16.5 










0/6,3 


51101 


17 


16/23.5 


0/2.3 












10/8.7 










0/2,8 


51145 


16 


0/3.0 


10/14.8 


16/12.6 


?/9.0 








16/18.1 






0/2.3 






51146 


16 


3/6.8 


0/4.9 




0/3.4 








2/6.6 












51147 


8 


0/3.5 




0/5.5 




















6/11,6 


51148 


5 


1/3.6 


5/31.2 












0/13.4 










0/6,4 


51149 


16 


0/5.4 


0/5.1 






0/2.4 




0/2.6 


13/11.4 


0/2.9 


0/2.1 




0/2.3 


0/2.5 1/5.2 


51150 


16 


0/5.4 




0/2.8 








0/2.2 


14/22.7 


0/2.3 






0/4.0 




51151 


9 


2/7.7 


9/15.1 












2/5.3 












51152 


17 


1/6.2 


4/7.6 




0/2.8 






0/2.5 


16/14.5 








0/4.6 


3/7,3 


51153 


16 


0/3.0 


0/4.0 


0/2.7 




0/5.4 


1/4.8 




15/16.1 






0/2.0 






51154 


19 


19/69.6 


5/9.0 








6/8.9 


0/2.3 


12/14.5 








0/2.7 




51155 


16 


14/18.1 


0/2.1 


0/4.8 






0/2.5 


0/2.9 


11/9.9 






0/6.2 


1/6.0 


0/2,2 


51156 


16 


10/49.0 


3/6.1 


4/14.4 


16/32.0 








16/52.2 












51157 


6 


0/6.9 


6/23,2 










0/3.2 


6/53.2 












51211 


16 


0/5.5 


2/6.4 


15/20.3 


4/6.4 


0/2.0 


0/4.6 


0/3.2 


0/2.2 


0/2.2 






1/4.1 


0/2,7 


51212 


16 


3/5.6 


0/3.8 




0/2.1 


0/2.9 




0/2.1 


16/38.6 






0/2.7 


0/2.3 


0/2.8 


51213 


16 




4/6.0 


0/2.6 


0/2.1 


0/7.1 




0/2.8 


0/2.8 








15/13.7 


13/14.1 


51214 


16 


2/5.9 


0/2.5 




0/2.9 




0/2.1 


0/3.7 


16/18.5 








0/3,3 




51215 


16 


3/6.2 


0/2.1 


0/2.5 


15/13.2 


0/2.1 


0/3.0 


0/2.1 


2/5.4 








11/12,4 




51216 


16 


9/9.7 


1/6.1 


10/10.0 


9/9.4 






0/2.8 


9/10.4 








15/19,7 




51217 


16 


0/4.0 


0/4.3 


0/3.6 




0/2.7 


7/8.4 


0/2.8 


14/28.8 








0/3,1 


0/2,9 


51218 


16 


2/8.6 


0/2.6 


0/3.9 


0/5.0 






0/2.8 


16/24.1 








0/4,7 


4/7,0 



Note, — For Hch observing dale we lisi ihe number of 17,5 min iniegraiions ihal were independenily seiirched. For each pulsar on each day we lisi ihe number of iniegraiions in which die pulsar was delected wiih Ihe search code according to ihe criteria 
summarized in §|t and Ihe si gnal-lo-noise ratio (cr) scaled to a 17,5 min integration, obtained from suimning the entire day's data set using the best pulsar ephemeris available (e.g., on MJD 51004, 47 Tuc C was not detected in any of 16 integrations with cr > 9, 
but all data addeff together using the pulsar's ephemeris yielded cr = 23.2, or an average of cr = 5.8 per integration. Hence the 0/5.8 notation. For 47 Tuc D, even adding all data did not yield an unambiguous detection — one with cr > 7 ■ — and so the respective 

entry is blank). 
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Table 3 

Pulse periods, dispersion measures, duty-cycles, and flux densities at 20 cm. 



Pulsar 


Period 


DM 


"'50 


Flux density 




(ms) 


(cm""^ pc) 


(%) 


(mJy) 


Previously Known 


47 Tuc C 


5.7567S0 


24.6 


13 


0.36(4) 


47 Tuc D 


5.357573 


24.7 


10 


0.22(3) 


47 Tuc E 


3.536329* 


24.3 


25 


0.21(3) 


47 Tuc F 


2.623579 


24.4 


14 


0.15(2) 


47 Tuc G 


4.040379 


24.4 


13 


0.05(2) 


47 Tuc H 


3.210341* 


24.4 


13 


0.09(2) 


47 Tuc 1 


3.484992* 


24.5 


14 


0.09(1) 


47 Tuc J 


2.100634* 


24.6 


15 


0.54(6) 


47 Tuc L 


4.346168 


24.4 


15 


0.04(1) 


47 Tuc M 


3.676643 


24.4 


25 


0.07(2) 


47 Tuc N 


3.053954 


24.6 


11 


0.03(1) 


Newly Discovered 


47 Tuc 


2.643343' 


24.4 


12 


0.10(1) 


47 Tuc P 


3.643021* 


24.3 


15 




47 Tuc Q 


4.033181* 


24.3 


11 


0.05(2) 


47 Tuc R 


3.480463* 


24.4 


15 




47 Tuc S 


2.830* 


24.4 


9 




47 Tuc T 


7.589* 


24.6 


17 




47 Tuc U 


4.342827* 


24.3 


8 


0.06(1) 


47 Tuc V 


4.810* 


24.1 


7 




47 Tuc W 


2.352344* 


24.3 


17 





Note. — Uncertainties in barycentric period and DM are one or less in the last digits quoted. 
*Binary pulsar (see Table a). 



the binary motion is to cause a change in the apparent pulse 
frequency during the integration, spreading the emitted signal 
power over a number of spectral bins, thereby reducing the ap- 
parent signal-to-noise ratio. In order to recover the signal, it is 
necessary to transform the time series to the rest frame of an 
inertial observer with respect to the pulsar before carrying out 
the periodicity search. This transformation is readily achieved 
by applying the Doppler formula to relate a time interval in the 
pulsar frame, r, to the corresponding interval in the observed 
frame, f : 

r(f) = To(l + v(f)/c), (1) 

where v{t) is the observed radial velocity of the pulsar along 
the line-of-sight, c is the speed of light, and we have neglected 
terms in (v/c) of order higher than the first. Ideally, if the orbital 
parameters of the binary system are known, v(f) can be calcu- 
lated from Kepler's laws. For the purposes of a blind search, 
where the orbital parameters are a priori unknown, assuming 
a Keplerian model for v(f) would require a five-dimensional 
search of all the parameter space — a computationally non- 
trivial task! The search can, however, be considerably sim- 
plified by assuming that the orbital acceleration a is constant 
during the observation, i.e., v(f) = at . This assumptio n is com- 
monly used in searches of this kind (Anderson 1992) and turns 
out to be reasonable over a wide range of orbital parameters, 
phases, a nd ob servation lengths. We investigate this in more 
detail in § |5.1.l| 

Given a prescription for v(f), the re-sampling process is 
straightforward. The time intervals in the new frame are cal- 
culated from equation |l} New sample values are then created 
based on a line ar interpolation running over t he original time 
series (see also Vliddleditch & Kristian 1984). We choose to 



pling interval, 
condition v(f ) : 



samp, at the mid-point of the integration. For the 
: at, it follows that 



fs; 



T0 = 



i+flr/2c' 



(2) 



define the value of tq such that r is equal to the original sam- 



where T is the integration time, 17.5 min in our case. This con- 
dition guarantees that, under the approximations made here, the 
number of samples in the corrected time series is the same as in 
the original one. The re-sampling process is relatively trivial in 
terms of computation time and can be readily carried out for a 
number of trial accelerations for each de-dispersed time series 
in order to effect a search in "acceleration space". 

In all the analyses carried out in this paper, each de-dispersed 
time series was re-sampled and analyzed for a number of trial 
accelerations in the range \a\ < 30m s~^. Most pulsars in 
presently known binary systems typically undergo accelera- 
tions of up to about 5 ms"^. There are however other binary 
systems where more extreme accelerations have motivated us 
to search a larger range of acceleration space. These include 
the double neuti'on-star binary PSR B2127H-I IC in M15, orig- 
inally discover ed at a trial acceleration of -9.5 ms~^ ( Ander- 
son et al. 1990), and the ec lipsing binary PSR B1744-24A in 
Terzan 5 dLyne etal. 1990| ), where the maximum line-of-sight 
acceleration is 33 m s"^. 

The step size between each trial acceleration was 0.3 ms"^. 
Some care is required when choosing this interval in order to 
strike a compromise between unnecessary processing time in- 
curred by oversampling the parameter space and the loss of sen- 
sitivity caused by undersampling. Assuming that the constant 
acceleration approximation adequately describes the Doppler 
shifts during the integration, it follows that the number of spec- 
tral bins across which a signal will drift due to the assumed trial 
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acceleration being different from the true value by an amount 
Afl is given by Aa x /Pc, where P is the pulse period. Our 
choice of step size thus guarantees that any pulsar period or 
harmonic with a period greater than 2 ms will not drift by more 
than one spectral bin during the 17.5 min integration. 

Following the de-dispersing and re-sampling procedures out- 
lined above, each time series is analyzed independently accord- 
ing to standard pulsar search techniques. The analysis proce- 
dure is very similar to that described in detail by Manchester 
et al. (1996). In brief, we compute discrete Fourier transforms 
of each time series, and sum harmonically related spectral com- 
ponents for 1, 2, 4, 8, and 16 harmonics in turn, in order to 
maintain sensitivity to a variety of pulse profile relative widths, 
or duty-cycles. For each harmonic sum, candidates are ranked 
in order of spectral signal-to-noise ratios^. At this stage, all 
spectral features with cr > 6 are stored, along with the period of 
the signal. 

Once computations are completed for the time series cor- 
responding to all trial accelerations for a 17.5 min block of 
data, we plot a versus trial acceleration for multiple occur- 
rences of each period, within a suitable tolerance, detected 
above a pre-defined "significance" threshold. This threshold 
can be estimated theoretically from the statistics of the phase- 
space searched, but practical considerations (such as the pres- 
ence of radio-frequency interference) result in our determining 
this value empirically. We generate plots such as that in Fig- 
ure |l] for all cases where at least one detection of a multiply 
detected period has cr > 9. 



detectable. Following Dewey et al. (1984), this is given by 
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1 
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-ZO 20 

Acceleration (m pr^) 

Fig. 1 . — Signal-to-noise ratio (cr) versus trial acceleration for the discovery 
of 47 Tuc R. This plot coiresponds to one of seven 17.5 min integrations, out of 
a total of eight on MJD 50742, in which a in at least one trial acceleration was 
higher than the threshold of 9. Most points with o" < 8 correspond to random 
noise rather than detection of the pulsar. 



y/nAiyT ^P 



1/2 



(3) 



Here cr is the signal-to-noise ratio threshold of the survey (9 
in our case); 77 is a constant that takes account of losses in 
sensitivity due to digitization of the incoming data and other 
inefficiencies in the system (~1.5 for our system); 5sys is the 
system equivalent flux density (^35 Jy for the present system 
in the direction of 47 Tuc); n is the number of polarizations 
summed (two in our case); Aiy is the total observing bandwidth 
(288 MHz); T is the integration time (1049 s); and w and P are 
the pulse width and period respectively. 

The observed pulse width w in the above expression is in 
general greater than the intrinsic width Wint emitted at the pulsar 
because of the dispersion and scattering of pulses by free elec- 
trons in the interstellar medium, and the post-detection integra- 
tion performed in hardware. Interstellar scattering is negligible 
for most observations toward 47 Tuc, which is relatively nearby 
and located well away from the Galactic plane. The observed 
pulse width for unaccelerated pulsars will therefore be the con- 
volution of the intrinsic pulse width and broadening functions 
due to dispersion and integration, and can be estimated from the 
quadrature sum 



DM ADM + 'samp- 



(4) 



Here foM is the dispersion broadening across one filter bank 
channel, Jadm is the overall pulse broadening due to the differ- 
ence between the single DM at which we de-disperse all data 
and the true DM of a pulsar (ADM), and fsamp is the data sam- 
pling interval (125 fis in our case). The dispersion broadening 
across a single channel is calculated using the standard expres- 
sion 



tuM - 8.3 



DM 

cm""* pc 



Ai^H 



MHz 



GHz 



MS, 



(5) 



which is valid for Ai/chan <C as in our case where Aj/chan = 
3 MHz and h' = 1374 MHz. The effective time resolution for 
these 47 Tuc observations, calculated from equation |] and the 
sample time, at a DM of 25 cm""* pc, is approximately 0.3 ms at 
the center frequency of observations, although it falls to 0.4 ms 
at the low-frequency end. However in many cases the effective 
time resolution is limited by the f adm term in equation H, which 



amounts to 0.9 ( 



ADM 



)msat 1374 MHz. 



Figure ^ shows the search sensitivity as a function of period 
for a 17.5 min integration based on equations ^|-^. Here we 
have assumed a box-car pulse shape and an intrinsic duty-cycle 
of 15%, typical of many among the pulsars known in 47 Tuc. It 
can be seen that the present system is sensitive to most pulsars 
in 47 Tuc with P > 2 ms and with average 20 cm flux densities 
above 0.3 mJy during a 17.5 min integration. 

^Spectral signal-to-noise ratios are defined throughout this paper to be the height of a feature in the amplitude spectrum divided by the spectral rms. In general, 



3. SENSITIVITY 

To estimate the sensitivity of our search we calculate the min- 
imum flux density Smin that a pulsar must have in order to be 



,latl 



dSxpm reconstructed time-domain pulse profiles are different, although in practice, for relatively low values 
^5.1.1) — they are identical, and we use the nomenclature a in both instances. To estimate average flux densities 



spectral signal-to-noise ratios and ratios call 
rplpygnt for Search sensitivity calculations ( 
(§p.fi), where some signal-to-noise ratios are Targe, we carry out all computations in the time domain. 




two of the 1 1 pulsars detected on at least 25% of the observing 
days. In contrast, at 50 cm, all but 47 Tuc C are detected less 
than 25% of th e time, and the same is true at 70 cm for seven 
of the pulsars (Robinson 1994). The detection rates with the 



Pi.ilse period (ms) 



Fig. 2. — Sensitivity as a function of pulse period. Tlie solid curve is 
the sensitivity calculated for our 20 cm observations, for an assumed pulsar 
duty-cycle of 15%. analyzed in individual 17.5 min blocks. The steps in the 
curves going from right to left reflect loss of power in the frequency domain 
as higher-order harmonics drop out of the spectrum. The dotted curve is an es- 
timate of the sensitivity at 50 cm for the observations reported by Manchester 
et al. (1990, 1991), for an integration time of 42min. The dashed curve scales 
the 20 cm curve to a wavelength of 50 cm assuming a pulsar spectral index of 
-1.6. 

The improvement in sensitivity provided by the new 20 cm 
system is apparent when comparing its sensitivity curve to that 
based on a 42 min observation assuming the parameters of the 
50 cm system used by Manchester et al. (1990, 1991) to dis- 
cover pulsars 47 Tuc C-M, also shown in Figure ^. For flat 
spectrum sources one may compare the 20 and 50 cm curves 
directly. A more realistic comparison between the two searches 
is shown by the scaled version of the present 20 cm curve as- 
suming the typical pulsar spectral index of -1.6. With this as- 
sumption we see that the present system represents as much as 
a threefold improvement in sensitivity for a 2 ms pulsar such 
as 47 Tuc J. 

The sensitivity calculations described in this section take no 
account of the effects of the degradation in a of any binary 
pulsar whose apparent period may change significantly during 
t he 17 .5 min integration. We investigate this issue in detail in 



§5.1.1 



4. RESULTS 

We have reduced all data listed in Table ^ in the manner 
described in §^. In this table we list the number of 17.5 min 
integrations analyzed on each day, giving both the number of 
daily integrations in which the search code detected the pulsar 
— i.e., in which (Tsearch > 9 — and the average daily cr, scaled 
to a 17.5 min integration, obtained from summing all data for 
that day using the pulsar's ephemeris. We detected some 50 
distinct periods for which the a vs. trial acceleration plots are 
broadly similar to that of Figure |l], i.e., have a clearly defined 
and significant maximum in trial acceleration. Eleven of these 
we identified with previously known pulsars in 47 Tuc. 

4. 1 . Previously Known Pulsars 

The detection rate for the previously known pulsars ranges 
from 10% (47 Tuc N), to 90% (47 Tuc C and J), with all but 



new observing system are clearly much higher than with past 
observations, leading to improved parameters for most pulsars. 

In Table ^we present "phase-connected" timing solutions for 
three previously known pulsars, obtained from 20 cm data col- 
lected since 1997 October To determine these we first make 
time-of-arrival (TOA) estimates by cross-correlating pulse pro- 
files, obtained by summing several minutes' worth of data 
for a pulsar on a given day, with a high-cr "standard pro- 
file", and adding the time delay to the start-time of the inte- 
gration, appropriately translated to the mid-point of the inte- 
gration. We then use the TOAs, and a model describing the 
pulsars, including astrometric, spin, and (where relevant) bi- 
nary parameters, to obtain improved parameters. These are de- 
rived by minimizing the timing residuals (difference between 
observed and computed pulse phases) in a least-squares sense. 
This procedure is implemented in the TEMPO softw are package 
(http://pulsar.princeton. edu/ tempo). 



The timing solutions for 47 Tuc C and D match those first 
presented by Robinson et al. (1995), obtained from 50 and 
70 cm data. The phase-coherent solution for 47 Tuc J is a new 
one. The small limit on its orbital eccentricity is perhaps some- 
what surprising, since it suggests a lack of strong recent in- 
teractions with other stars, while its position well outside the 
cluster core tends to suggest the opposite. However it is also 
conceivable that tidal dissipation in a Roche lobe-filling com- 
panion could circularize the orbit quickly (cf. §5.4). We remark 
in passing that a 20 cm continuum map of the central regions 
of 47 Tuc shows three re latively bright, variable, point sources 
(Fruchter & Goss 1999). Two of these are clearly associated 



with 47 Tuc C and D. The third, whose position differs from 
that Hsted for 47 Tuc J in Table | by 0?07(3) in right ascen- 
sion and 0."3(2) in declination, and whose average flux density 
of 0.6 mJy matches our average of 0.5 mJy (§5.6), can now be 
identified with 47 Tuc J. 

There is no evidence that 47 Tuc L is a binary pulsar (a pos- 
sibility according to Robinson et al. 1995). We have detected 
this pulsar on 1 1 separate occasions spanning over 500 days 
(Table and the barycentric period is constant (see Table ^ 
which lists the pulse periods, DMs, and observed duty-cycles 
of all pulsars detected in our observations). 

Finally, we were also able to determine the eccentricity for 
47 Tuc E, and obtain a binary solution for 47 Tuc H. The very 
significant eccentricities for these systems (Table |]) are unique 
among this class of b inaries: in the disk of the Galaxy such sys- 
tems have e <; lO"*" ( [Camilo 1999 ). The eccentricities are pre- 
sumably fossil evidence of interactions near the dense cluster 
core, and may be due to pert urbations by passing st ars, or stel- 
lar exchange i nteractions (see|Rasio & Heggie 199^ and H eggie 
& Rasio 1996 for details). It will be interesting to see if other 
wide systems in this cluster have large eccentricities. 

4.2. Newly Discovered Pulsars 

Nine of the remaining periods detected in our data appeared 
in at least two 17.5 min integrations on one day, with some be- 
ing detected in most integrations. For each detection of such 
a period, we de-dispersed the raw data at several trial values 
of DM about the nominal DM. In all nine cases there was a 
clearly defined maximum signal-to-noise ratio as a function of 
DM, ensuring that the periodic signals were truly dispersed. Af- 
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Table 4 

Timing solutions of three pulsars. 



Parameter 


J0023-7204C 


J0024-7204D 


J0023-7203J 


Right ascension, a (J2000) 


00''23'"50?350(1) 
-72°04'31"486(3) 


00''24"'13'877(2) 
-72°04'43"845(9) 


00''23"'59!4057(6) 
-72°03'58"781(2) 


Declination, 5 (J2000) 


Period, P (ms) 


5.7567800118(2) 


5.3575732858(1) 


2.100633545861(1) 


Period derivative, P 


-5.1(1) X 10"^° 


-0.31(4) X 10"^° 


-0.96(1) X 10"-° 


Epocii of period (MJD) 


47858.5 


48040.7 


51000.0 


Dispersion measure, DM (cm""' pc) . . . 


24.6 


24.7 


24.6 


Orbital period, Pj, (d) 






0.120664939(1) 


Projected semi-major axis, .v (s) 






0.040411(1) 


Eccentricity, e 






< 0.0001 


Longitude of periastron, u 






0.0 


Time of periastron. To (MJD) 






51000.765678(1) 



Post-fit rms timing residual (/is) 
a - Q47 T,ic 

(5 - 547 Tuc 

Pulsar-cluster center offset (core radii) 



10 

-15?6 
19"6 
6 



19 

8!0 
7"3 
3 



9 
-6^5 
52"3 
5 



Note. — Figures in parentheses represent uncertainties in the last digits quoted, given for all fitted parameters. 



ter viewing the corresponding time-domain pulse profiles, and 
judging them suitably "pulsar-like", we felt secure in identify- 
ing each of these nine periods with previously unknown pulsars 
in 47 Tuc. Table ^ lists detections for the newly discovered 
pulsars that are seen only occasionally. 

For all nine of the newly detected pulsars, the period changed 
within an observing day in a manner consistent with the 
Doppler shifts arising in binary systems. Only three of the nine 
pulsars (47 Tuc O, Q and U) were detected on many occasions 
(Table g). Fortunately, three of the remaining six (47 Tuc P, 
R, and W) have such short orbital periods that it was possible 
to cover more than an entire orbit on their discovery days, and 
hence establish good orbital parameters (see Fig.Jsl). Binary 
parameters for all pulsars are summarized in Table ^ 

Most of the newly discovered pulsars appear to be similar 
in their binary characteristics to the previously known pulsars 
in the cluster. For instance, 47 Tuc O, P, and R have short 
orbital periods (Pt, ^ 1.5-5.5 h) and low companion masses 
(m2 ^ 0.03 Mq), similar to the previously known 47 Tuc I and 
J. Although the period of 47 Tuc R, at 96min, is considerably 
shorter than that of 47 Tuc J, there is no reason to suspect that 
these systems are fundamentally different. 

In contrast, 47 Tuc Q and U are similar to 47 Tuc E and H, 
with OT2 ~ 0.2 M0; the companion stars are presumably low- 
mass helium white dwarfs much like those found in similar 
syste ms in the Galactic disk (e.g., L undgren, Foster, & Camilo 
1996). The binary periods for 47 Tuc Q and U are somewhat 
shorter than the corresponding ones for 47 Tuc E and H (0.4- 
1 .2 d vs. 2.3 d; see Table ||), but we believe this betrays no fun- 
damental difference between the various systems. We have not 
yet obtained orbital parameters for 47 Tuc S and T, but from the 
observed projected velocity change during the discovery days 
(-5 and 12kms~' over 2 h, respectively; see Fig. ||) we believe 
they also fall into this higher-mass group. 

It thus appears the vast majority of the binary pulsars known 
to date in 47 Tuc fall into one of two classes. However, the two 
pulsars discovered most recently, 47 Tuc V and W, are not so 
easily categorized. 




Time (h) 

Fig. 3. — Ob.served barycentric periods for six newly discovered pulsars 
that are seldom detected. Each period is based on 17.5 min of data, with the 
exception of 47 Tuc R (8.7 min) and 47 Tuc V (2.2 min), arid corresponds to 
<j >6. Times are given in UT on discovery day (see Table H). 

Pulsar 47 Tuc V has been detected on two different days, at 
seemingly identical orbital phases (with period increasing as a 
function of time, hence on the "near-side" of its companion — 
see Fig. ||), and it displays an apparent orbital velocity change 
of up to ~ 80km s"' in under one hour The observed change 
in period during this time leads us to believe that the orbital 
period is several hours, suggesting a relatively massive com- 
panion. Another unusual aspect of this system is that it was 
detected on both occasions with relatively large signal-to-noise 
ratios in several 2 min sub-integrations, but interspersed with 
non-detections (Fig. as if it were being irregularly eclipsed, 
despite its orbital phase^. 

''in practice it is usually easy to distinguish among interstellar scintillation and eclipses by a companion as the cause for the temporary non-detection of a pulsar 
Eclipses tend to occur with the pulsar roughly centered at superior conjunction, while scintillation has no preferred orbital phase. Also, scintillation has characteristic 
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Table 5 
Binary parameters. 



Pulsar 




.V 


To 


u 


e 






(d) 


(s) 


(MJD) 


(deg) 




(Mq) 


Previously Known 


47 Tuc E 


2.2568448 


1.98184 


51004.05 


219 


0.0003 


0.19 


47 Tuc H 


2.3576966 


2.15281 


51000.9735 


110 


0.071 


0.20 


47 Tuc I 


0.2297923 


0.03846 


50740.5794 


0.0 


0.0 


0.02 


47 Tuc J 


0.1206649 


0.04041 


51000.7657 


0.0 


0.0 


0.03 


Newly Discovered 


47 Tuc O 


0.1359743 


0.04515 


51000.0212 


0.0 


0.0 


0.03 


47 Tuc P 


0.1472 


0.0380 


50689.6789 


0.0 


0.0 


0.02 


47 Tuc Q 


1.1890840 


1.46219 


51002.175 


0.0 


0.0 


0.21 


47 Tuc R 


0.0662 


0.0334 


50742.6365 


0.0 


0.0 


0.03 


47 Tuc S 


t 












47 Tuc T 


t 












47 Tuc U 


0.4291057 


0.52696 


51002.645 


0.0 


0.0 


0.15 


47 Tuc V 


~0.2? 












47 Tuc W 


0.1330 


0.2435 


51214.950 


0.0 


0.0 


0.15 



Note. — Uncertainties in all fitted parameters are one or less in the last digits quoted. 
* Assumes a pulsar mass of 1 .4 Mq and inclination angle for the orbit of / = 60° . 
^ About one to several days. 



Pulsar 47 Tuc W, with an orbital period of 3.2 h, may resem- 
ble short-period binary systems such as 47 Tuc J, and the pulsed 
radio signals are apparently eclipse d fo r a fraction of its orbit, 
as in some of these systems (see §5.4). However its implied 
companion mass, at ~ 0.15 M0, is a factor of ^ 6 higher than 
that of those systems, more typical of a low-mass white-dwarf. 
But if the companion were indeed a white dwarf, we would 
not expect eclipses. We consider it possible that this system 
is somewhat akin to PSR B1744-24A, located in the cluster 
Terzan 5, with Pf, ^ 1.8 h, OT2 ^ 0.10M(^, and display ing irreg- 



ular eclipses (Lyne et al. 1990; Mice & Thorsett 1992) 



There remain about 30 periodic signals that displayed sev- 
eral pulsar-like characteristics, but that were detected on only 
one integration of one day. We regard several of these as good 
pulsar candidates, potentially to be confirmed by future obser- 
vations. 

5. DISCUSSION 

5.1. Survey Limitations 

Although the recent observations of 47 Tuc have been suc- 
cessful beyond our most optimistic expectations, it is important 
to keep in mind the many and severe selection effects inherent 
in the data and analyses reported on in this paper. The greatest 
of these, particularly for short-period binary systems, has to do 
with the extent of parameter space that has been covered by our 
"acceleration searches". 

5.1.1. Sensitivity to Binary Pulsars 

In addition to the improvements in sensitivity over previous 
searches of 47 Tuc enabled by the new 20 cm observing system, 
use of an acceleration search has significantly improved the 

time-scales and amplitudes. No other pulsar in 47 Tuc displays the deep intensity modulations on time-scales of about a minute observed for 47 Tuc V, and we 
therefore suggest the eclipse interpretation. 

'in effect this suggests that without the acceleration code we would have discovered five new pulsars, an increase of ~ 50% over the previously known population. 
This is surprising, considering the factor of ~ 3 improvement in sensitivity claimed in §bI. This estimate relied on assuming a spectral index of —1.6 for pulsars in 
47 Tuc. Perhaps this assumption is not valid, and the average spectral index for millisecotid pulsars in 47 Tuc is smaller than —1.6. 



ability to find binary pulsars. Prior to this survey only four out 
of 1 1 (36%) millisecond pulsars known in 47 Tuc were mem- 
bers of binary systems. The new discoveries bring the number 
of binary pulsars in 47 Tuc to 13 (65% of the total). We con- 
sider it highly significant that all of the new millisecond pulsars 
reported here are members of binary systems. These discov- 
eries are in spite of the fact that, as demonstrated in §^ the 
sensitivity to isolated pulsars is significantly improved over the 
previous surveys by Manchester et al. (1990, 1991) and Robin- 
son et al. (1995). We note that three of the five brightest pulsars 
in 47 Tuc at 20 cm are isolated (Table ||); if significant, this may 
suggest that isolated pulsars in 47 Tuc have a relatively high in- 
trinsic luminosity cutoff. Also, isolated pulsars in 47 Tuc might 
have, on average, a larger spectral index than binary pulsars. Fi- 
nally, the specific incidence of isolated pulsars in 47 Tuc might 
be low. A combination of these factors could explain why our 
more sensitive survey has failed to find additional isolated pul- 
sars. 

The improved sensitivity of the acceleration search to short- 
period binaries is demonstrated by the fact that four of the 
nine binary systems we found, 47 Tuc P, R, S and V, would 
not have been discovered without the use of the acceleration 
code since, in each case, the signal-to-noise ratio never rose 
above the threshold of 9 at a trial acceleration of zero^. Despite 
searching an acceleration range of |a| < 30 ms"-, most of the 
new discoveries and detections of previously known pulsars oc- 
curred at |fl| < 5 ms~^. Three notable exceptions are 47 Tuc R, 
V, and W, that displayed maximum line-of-sight accelerations 
of 11.4, 23.7, and 21.3ms"^ respectively (e.g.. Fig. |l]). 

While it is clear that this survey has improved the knowledge 
of the true fraction of binary pulsars in 47 Tuc, an outstanding 
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Table 6 
Rarely detected pulsars. 



Pulsar 


Date 


Number of 


Detections/o" 




(MJD) 


integrations 




47 Tuc P 


50689 


16 


5/6.3 


47 Tuc R 


50742 


8 


7/8.9 


47 Tuc S 


50741 


8 


3/9.9 


47 Tuc T 


50746 


8 


4/9.3 




51005 


8 


0/5.4 




51040 


4 


4/25.8 


47 Tuc V 


51012 


8 


2/3.4 




51055 


16 


3/4.4 


47 Tuc W 


51214 


16 


9/10.0 



question that remains is the extent of binary parameter space 
probed by the present search strategy. In this section we attempt 
to answer this by exploring the limits of the present search with 
respect to a variety of short orbital period binary systems. In 
what follows, when we use the term "acceleration search" we 
are explicitly referring to the one-dimensional (i.e., constant ac- 
celeration) approximation outlined in §^ 

For the purposes of this discussion it is useful to note that 
most observed binary pulsar systems in 47 Tuc fall broadly into 
two categories: those with orbital periods of order 0.4-2.3 d and 
companion masses ~ 0.2 M© (47 Tuc E, H, Q, U, and probably 
S and T) and 47 Tuc I, J, O, P, and R, which are characterized 
by shorter orbital periods (1.5-5.5h) and lighter companions 
(m2 ~ 0.03 Mq). (The nature of 47 Tuc V and W is unclear at 
present.) We shall refer to these types here respectively as the 
"normal" and "short-period" binary systems, investigating the 
sensitivity of our search to both types in turn. In addition, since 
all the presently known binary systems have essentially circu- 
lar orbits, we also investigate the search sensitivity to putative 
binary systems in which the pulsar is in an eccentric orbit about 
a more massive companion. 

The first, and really only, detailed study of the degradation 
in sensitivity of a radio pulsar search for short orbital period 
systems was carried out by Johnston & Kulkarni (1991; here- 
after JK). By calculating the decrease in amplitude of the signal 
in the Fourier power spectrum for a binary pulsar relative to a 
solitary one, and mathematically describing the signal recovery 
technique involved in the acceleration search, JK were able to 
quantify the degradation in sensitivity of pulsar search codes 
with and without acceleration schemes of the type used here. 
As their main example, JK considered a 1000 s observation 
of a 1.4M0 pulsar in a circular orbit with a 0.3 compan- 
ion for a wide variety of orbital and pulse periods. Given that 
our 17.5 min analyses correspond to a similar integration time, 
JK's results are highly relevant to this discussion. In addition, 
the 0.3 M0 companion considered by JK closely resembles the 
class of normal binary systems defined above. 

We focus on Figure 4 of JK, which shows the net gain in 
sensitivity of an acceleration search when making multiple ob- 
servations of the same system. In their example, JK assumed 
five independent observations of the binary at random orbital 
phases. From this figure we see that, for a 10 ms pulsar, the ac- 
celeration search allows one to detect binary systems with or- 
bital periods of order four times shorter than a standard pulsar 
search code. The gain increases still further for shorter period 
pulsars, where the blurring of harmonics in the power spectrum 



becomes even more pronounced. For a 3.5 ms pulsar, similar to 
47 Tuc E and H, we infer from Figure 4 of JK that the accel- 
eration analysis allows the detection of a binary with an orbital 
period about six times shorter than the standard search. 

The results of JK are not directly applicable to the popula- 
tion of short-period binaries in 47 Tuc. Presently, the 3.48 ms 
pulsar 47 Tuc R, with an orbital period of 96 min, is the shortest 
orbital period system of this type. In order to estimate the sen- 
sitivity of our search code to even shorter period binaries with 
low-mass companions like 47 Tuc R, we generated a number of 
time series containing a synthetic 3.48 ms pulsar whose period 
was modulated as in a circular Keplerian orbit with a low-mass 
companion. We considered four separate cases corresponding 
to orbital periods of 90, 60, 30 and 15 min. The first value con- 
sidered is similar to that of 47 Tuc R; the last value is close 
to that of the shortest orbital period known for any binary sys- 
tem: the 1 1 min orbit of the X-ray burst source XI 820-303 in 
the g lobular cluster NGC 6624 (S tella, Priedhorsky, & White 
1987). 

In each case the mass function was held constant at 8.68 x 
10"^ M0 (similar to that of 47 Tuc R). To quantify the im- 
provement gained by using the acceleration code, we generated 
a control time series containing a fake 3.48 ms pulsar with a 
fixed period but otherwise identical to the other simulations. In 
a similar sense to JK, we define the efficiency factor 7 as the 
detected a of the binary pulsar divided by that of the control 
pulsar In order to compare the relative efficiency of the non- 
accelerated and accelerated analyses, separate values of 7 were 
computed: 70 is the efficiency for the non-accelerated analysis, 
and jA is the equivalent parameter for the acceleration search 
in the analysis of the same data, performed as described in §^ 
By definition < 7 < 1, with a value of unity signifying that 
no loss of signal occurred. In each case, the simulations were 
performed over all orbital phases to establish the response of 
the search code to the full range of orbital accelerations. 

The results of the simulations are summarized in Table ^ 
where we list the range of efficiencies, as well as the mean 
values averaged over all orbital phases. As expected, the ef- 
ficiencies for both the normal (70) and acceleration search code 
(7a) are strong functions of orbital phase. In Table ^ we also 
list the "improvement factor" / = 7^/70, the net gain in sig- 
nal achieved by employing the acceleration search, averaged 
by time spent at each orbital phase. These simulations show 
a significant improvement in sensitivity gained when using an 
acceleration search. In particular for orbital periods as short 
as one hour (cases 1 and 2) the efficiency of the acceleration 
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Table 7 

Search code tests on simulated data. 



Case 


Period 


Pi, 


/ 




70.mi„ 


70 


To, max 


7/l,min 


7,1 




/ 




(ms) 


(min) 


(Mq) 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 




1 


3.48 


90 


8.68x10"" 


0.0 


33 


48 


75 


80 


88 


99 


1.8 


2 


3.48 


60 


8.68x10"" 


0.0 


26 


39 


58 


68 


75 


86 


1.9 


3 


3.48 


30 


8.68 xur" 


0.0 


24 


33 


41 


38 


45 


51 


1.4 


4 


3.48 


15 


8.68 xur" 


0.0 


18 


23 


28 


30 


32 


35 


1.4 


5 


34.8 


60 


0.13 


0.7 





48 


78 


18 


63 


92 


1.3 


6 


34.8 


30 


0.13 


0.7 





33 


63 


20 


43 


69 


1.3 



Note. — For each case considered, the columns list the assumed pulse period, orbital period, mass function, and orbital eccentricity. The remaining columns show the range and average 
values of efficiences of the search code for a normal analysis, 70, and using an acceleration search, 7^. The averaged "improvement factor" / = -^aIio is also given. Averaging over orbital 
phases has been appropriately weighted by the time spent at each phase. 



search is at least 70% for any orbital phase, more than twice the 
efficiency of the normal search code. 

For the shortest orbital periods of this type considered (cases 
3 and 4), the acceleration search continues to be more effi- 
cient than the normal search although the improvement factor 
is smaller. In these cases the integration is becoming an in- 
creasingly significant fraction of the orbital period. As a re- 
sult, the constant acceleration approximation used in the search 
(§||) ceases to hold. Among possible improvements to this 
strategy, the most straightforward approach is to perform one- 
dimensional acceleration searches on shorter segments of data. 
One can also include the acceleration derivative term in the pre- 
scription for v(t). Although this necessarily adds a dimension 
to the search, it may provide a dramatic improvement in cov- 
erage of orbital parameter space for radio pulsars with binary 
characteristics akin to XI 820-303. 

Finally, in Table we also summarize the results of simu- 
lations carried out to test the efficiency of our code at detect- 
ing a putative relativistic binary double-neutron star system in a 
highly eccentric short-period orbit (cases 5 and 6). As would be 
expected, for certain orbital phases a normal search has no hope 
of detecting such systems, while the acceleration search main- 
tains a detection efficiency of at least 20% at all orbital phases. 
Our search would therefore have found a relatively bright rel- 
ativistic binary pulsar in 47 Tuc at any orbital phase, but none 
was detected. This provides some empirical evidence that the 
population of pulsars in the central regions of 47 Tuc is com- 
posed exclusively of bona fide millisecond pulsars: those with 
rotation periods under ^ 10 ms, and in relatively circular orbits 
with low-mass companions, or with no companions at all. 

5.1.2. Other Selection Effects 

An interesting issue from the point of view of the origin and 
evolution of the pulsars in 47 Tuc is their intrinsic pulse pe- 
riod distribution. The observed distribution shown in Figure Q 
spans a remarkably narrow range and suggests that there may 
be a deficit of pulsars below 2 ms. This interpretation should 
be viewed with some caution, however, given the small number 
statistics, and it is clear from Figure ^ that we begin to lose sen- 
sitivity to pulsars with P <,2 ms. Presently, therefore, the ob- 
servations do not unambiguously address the question of what 
the shortest pulsar periods might be. Pulsar 47 Tuc W, for ex- 
ample, was originally detected at a period of 1.18 ms, showing 
that we are sensitive to some rather short periods, but to be truly 



sensitive to f <; 1 ms, a different experiment with higher time 
and frequency resolution will have to be carried out. 




4 6 
Pulse period (ms) 

Fig. 4. — The observed period distribution of the 20 pulsars known in 
47 Tuc. It is presently unclear to what extent the apparent deficit of pulsars 
wilh_oeriods below 2 ms is a result of period-dependent selection effects (see 
55.1.2). 



It is likely that the vast majority of all pulsars in 47 Tuc are 
to be found within the area covered by the 14'-wide 20 cm tele- 
scope beam pattern. This is in keeping with expectations based 
on the radial density profile of 47 Tuc ( |De Marchi et al. 1996 ) 
from which we infer that almost 70% of the stars within the 
tidal radius of the cluster (Table |l]) fall within the telescope 
beam area. This is additionally suggested by the observation 
that the pulsars for which a timing solution has been obtained 
are all located within a region with linear size ~ 2', and that 
we have re-detected all pulsars discovered at lower frequencies, 
even though the area subtended on the sky by the 20 cm beam 
pattern is at most one-sixth that within which the first 1 1 pul- 
sars were found. Nevertheless it is possible that some pulsars 
have been ejected from the cent ral regions of the cluster, as i s 
thought to be the case for M15 (Phinney & Sigurdsson 1991 ), 
and efforts to widen the search region might prove instructive. 
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There are a number of additional selection effects having to 
do with the manner in which we have reduced the data so far. 

For isolated pulsars, the sensitivity limits might be improved 
by a factor of up to ^ \/T6 by reducing the entire 4.66 h data 
sets in a coherent manner This could be done at reasonable 
computing expense because only the trivial zero-acceleration 
case would have to be addressed (once the Earth's rotation is ac- 
counted for). However, since flux densities for pulsars in 47 Tuc 
certainly vary significantly due to scintillation on time-scales of 
hours or less, the gain in sensitivity might be quenched by peri- 
ods within the 4.66 h during which we would be adding mostly 
noise to the time series being analyzed. In extreme cases, the 
sensitivity could actually be degraded! Optimally then, each 
session's data stream should be analyzed with a multiplicity of 
integration spans. This can, however, degenerate into a huge 
overall data reduction task. These considerations also apply to 
the case of binary pulsars, where the additional need to per- 
form sensitive acceleration searches over longer time intervals 
further exacerbates the overall computing load. 

Finally, as already noted, we have analyzed most of the data 
at only one value of DM. The DMs for known pulsars in 47 Tuc 
span about one unit (Table ||), and this akeady introduces some 
loss of signal for the shortest period pulsars we have consid- 
ered (§||). It would also be advantageous to broaden the range 
of DM searches because greater variance in the value of mea- 
sured DMs might add to the understanding of the properties of 
the ionized intra-cluster, and interstellar, medium. Most of the 
variance observed for DMs of pulsars in globular clusters is ex- 
pected to be due to slightly different lines-of-sight across the 
Gala xy, because m ost clusters are believed to be evacuated of 
gas (ISpergel 199 1|). 



5.2. Isolated Pulsars in 47 Tuc 

Mindful of the comments in the previous section, we can still 
make some general statements about the likely intrinsic popu- 
lation of pulsars in 47 Tuc. The observed population of pul- 
sars, all with millisecond periods, divides at present roughly 
into one-third isolated pulsars, one-third normal binary systems 
(Pt ^ 0.4-2.3 d and m2 ^ 0.2 M©), and one-third short-period 
binary systems (Ph ^ 0.06-0.2 d and m2 ^ 0.03 M©). 

As mentioned previously, the fraction of isolated pulsars 
known in 47 Tuc prior to this survey was 64% of the total. 
Use of acceleration code has improved the detectability of bi- 
nary systems, but we are still selecting severely against binaries 
of the short-period variety, when compared to isolated pulsars, 
and we consider it probable that the actual fraction of isolated 
pulsars in 47 Tuc is less than the 35% found at present. By 
comparison, the fraction of isolate d millisecond pulsars in the 
disk of the Galaxy is ^ 20% (e.g., Camilo 1999). The surveys 
of the Galactic disk, because of their short integration times, 
are relatively more sensitive to short-period binaries, so the iso- 
lated fraction for disk pulsars, up to small-number statistics, is 
already known with some confidence. 

The overall small fraction of isolated millisecond pulsars 
in 47 Tuc suggests an interesting puzzle: long-period bina- 
ries (Pt, ~ lO-lO-'d) are expected to be formed at high rates 
in dense clusters through exchange interactions between old 
isolated neutron stars and "hard" primordial binaries ( Hut et 
al. 1992; |Sigurdsson & Phinney 1993| ). The long-period bina- 
ries have, however, high cross-section to disruption, and should 
eventually leave behind isolated millisecond pulsars. We de- 
tect no long-period binaries, and relatively few isolated pulsars. 
Perhaps most wide binaiies ai^e disrupted before the neutron star 



companions can evolve into red giants, followed by mass trans- 
fer an d birth of milli second pulsars, as in the standard scenario 
(e.g., Verbunt 1993). This is an open question. In any case, 
even if the isolated fraction of millisecond pulsars in 47 Tuc and 
in the Galactic disk are comparable it is surely a coincidence, 
given the different evolutionary pressures present in both envi- 
ronments. 

5.3. Profile Morphology 

Chen & Ruderman (1993) suggested that the formation 
mechanism for millisecond pulsars in the Galactic disk and 
in globular clusters might be different. The observational ev- 
idence to support this idea relied on supposed differences in 
the observed radio beams (and hence magnetic field geome- 
try, and origin) when comparing the collection of pulse pro- 
files then available for disk millisecond pulsars and those in 
47 Tuc, and in particular the supposed dearth of interpulses 
among the pulsars in 47 Tuc. More recently, Jayawardhana & 
Grindlay (1996) also presented evidence suggesting that the ob- 
served pulse width distribution of disk and cluster millisecond 
pulsars may be different. 

Our collection of profiles for pulsars in 47 Tuc, shown in 
Figure ^, is larger and of higher quality than that presented by 
Manchester et al. (1991), on which Chen & Ruderman (1993) 
and Jayawardhana & Grindlay (1996) based their conclusions. 
In addition, there are now many more hig h-quality pulse pro 



files available for disk millisecond pulsars (Kramer et al. 1998 
Stairs, Thorsett, & Camilo 1999). Given these observational 



improvements, it is appropriate to briefly review the situation 
regarding possible morphological differences. In an attempt to 
ensure homogeneity, we choose to compare our 20 cm profiles 
with the sample of 19 Galactic disk pulsars reported by Kramer 
et al. ( 1998) which were also obtained from 20 cm observations. 
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Fig. 5. — Average pulse profiles at 20cm for pulsars in 47 Tu- 
canae as a function of rotational phase. These profiles are avail- 
able in digital form at the European Pulsar Network data archive 
(http://www.mpifr-bonn.mpg.de/div/pulsar/data). The ef- 
fective time resolution for all profiles displayed is ~ 0.3 ms, ranging from 0.04 
in pulse ph^e for 47 Tuc T to 0. 14 for 47 Tuc J. Some profiles are not resolved 
(see Table g). 

In order to test the suggestion that the pulse widths of the pul- 
sars in 47 Tuc are broader than their counterparts in the Galac- 
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tic disk, we compared the pulse width distribution measured at 
the 50% intensity level (wgo) for the 20 pulsars listed in Table |^ 
with that of the 19 disk millisecond pulsars presented in Table 2 
of Kramer et al. (1998). A Kolmogorov-Smirnoff test returns a 
25% probability that the two samples are drawn from the same 
underlying distribution. At best this is a marginally significant 
result, reflecting the larger fraction of Galactic disk millisec- 
ond pulsars for which w^q < 10% (but note that in spite of the 
improved time resolution of the new data, not all the profiles 
shown in Figure |]are fully resolved). 

Finally, from Figure |] we note that 1 1 pulsars have essen- 
tially simple, single-component profiles; five pulsars have clear 
interpulses; and four more have profiles that may be described 
as complex, with emission evident across much of the period or 
in several distinct components. Inspecting the sample of Galac- 
tic disk millisecond pulsars presented by Kramer et al. (1998), 
we find five clear interpulses from a sample of 24 pulsars. 
Given the present statistics, this is entirely consistent with the 
fraction of interpulses we observe in 47 Tuc, refuting the mor- 
phological differences suggested by Chen & Ruderman (1993). 

5.4. Eclipsing Pulsars 

Robinson et al. (1995) show that, at 70 and 50cm, 47 Tuc J 
is eclipsed by its companion for about a quarter of its orbit. We 
find that at 20 cm, on the other hand, the pulsar is always visi- 
ble at all orbital phases. However 47 Tuc O, a pulsar system in 
many respects similar to 47 Tuc J, is always eclipsed for about 
10-15% of its orbit at 20 cm. 

There are three other similar systems known in 47 Tuc, with 
orbital periods of less than 6 h and with ^ 0.03 Mq compan- 
ions. 47 Tuc P and R were detected at only one epoch, al- 
though at least one full orbit was covered in each case: 47 Tuc P 
was seen at all orbital phases, while 47 Tuc R was apparently 
echpsed for ^ 25% of its orbit, centered with the pulsar at supe- 
rior conjunction (Fig. ||). 47 Tuc 1 is detected much more often, 
and there are no eclipses observed at 20 cm. 

Finally, 47 Tuc V also displays apparent eclipses, but we 
cannot yet cha ract erize the system or its eclipses, which may 
be irregular (§ [4.2[ ); and 47 Tuc W, with m2 ~ 0.15 Mq, was 
eclipsed for ~ 25% of its orbit on the only occasion it was de- 
tected (Fig. 1^). While based on orbital parameters we suggest 
that the system may be somewhat similar to the eclipsing pulsar 
B1744-24A in Terzan 5 (§4-. 2), we cannot yet be sure. 

In the Galactic disk there are two pulsar systems known that 
have properties broadly similar to some of the above eclips- 
ing pulsars. PSRs B195 74-20 and J2051-0827 (A rzoumanian, 
Fruchter, & Taylor 1994; |Stappers et al. 1996] ) are millisec- 
ond pulsars in very circular orbits with ^ 0.03 Mq companions, 
with orbital periods of a few hours. These pulsars are both 
eclipsed by their companions, apparently bloated to fill their 
Roche lobes, but the systems differ in eclipse phenomenology: 
PSR B1957H-20 always displays eclipses for ~ 20% of its orbit, 
at all frequencies, while PSR J205 1-0827 is eclipsed at low ra- 
dio frequencies, but is visible at all orbital phases at 20 cm; in 
the latter instance, however, the radio pulses suffer extra time 
delays as they propagate in a dispersive medium. The differ- 
ences observed in these systems are likely due to a combination 
of different masses and chemical compositions of the compan- 
ions, and possibly evolutionary states; the energy that is de- 
posited at the companion by the relativistic pulsar wind, which 
in turn depends on the pulsar spin parameters, beaming geom- 
etry, and orbital separation; and geometrical effects. It is not 
yet clear if the echpse behavior for 47 Tuc J or O, let alone that 



for the seldom-detected pulsars in 47 Tuc, is similar to that of 
either PSR J2051-0827 or B1957H-20. 

5.5. The Central Mass Density of 47 Tuc 

In Table Q we list phase-connected timing solutions for three 
pulsars. An unusual feature of all the solutions is that the period 
derivatives, P, are negative, apparently implying the pulsars are 
spinning up, rather than spinning down as is the case for all 
pulsars known in the Galactic disk, which are losing rotational 
energy emitted in the form of magnetic dipole radiation and a 
relativistic particle wind. In fact this is a feature observed in 
som e other globular cluster pulsars (see, e.g., W olszczan et al. 
1989), and is most easily interpreted by considering that the 
observed Pobs is a combination of the intrinsic Pint(> 0) and a 
contribution resulting from the acceleration of the pulsar in the 
gravitational potential of the cluster, such that 



(6) 



In this equation, Pacc/P = a/fc, where a/ is the acceleration 
along the line of sight (Phinney 1992). A negative Pobs implies 
the pulsar lies behind the cluster center as seen from the Sun. 

Given a model for the mass distribution in the cluster, a lower 
limit on «/ by assuming Pint = 0, and an accurate pulsar position, 
one can estimate a lower limit to the central mass density in the 
cluster. Assuming that the central regions of the cluster can be 
described by King's (1962) empirical density model, one ob- 
tains for the central central mass density 
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In this expression r is the projected distance of the pulsar from 
the cluster center and k is the Jeans length at the cluster center. 
For a King model, we have made the standard approximation 
that K is one third of the core radius of the cluster, r^. . We use 
equation ^ for the pulsars in 47 Tuc by taking the most recent 
and reliable determination of the angular core radius, based on 
a fit of a King radial density profile to the observed stellar pop- 
ulation; together with a distance estimate to 47 Tuc (Table [ij), 
equation ^becomes 



p{Q) > 2.8 X 10^ 



arcmm 



10- 



-20 



ms 



Mqpc 



(8) 

Here R is the angular separation between the pulsar and the cen- 
ter of the cluster (see Table |l]). Using the celestial coordinates, 
P, and P, listed in Table || for 47 Tuc C, D, and J, we obtain 
from equation H lower limits on p{Q) of, respectively, 30, 1, and 
13 X IO^MqPC-I 

The most stringent lower limit on p{Q), which we obtain from 
the position and period derivative of 47 Tuc C, is some four 
times larger than the value quoted by Robinson et al. (1995) for 
the same pulsar. This is due to ou r use of the latest core radius 
estimate ( pe Marchi et al. 1996 ) which is approximately hal f 
the value of previous estimates (see, e.g., [Calzetti et al. 1993 ). 
De Marchi et al. demonstrate that the difference in core radius 
estimates can be entirely attributed to the effects of photometric 
incompleteness in the earlier samples. 

If a King model can be used to describe the central regions of 
the cluster, and the observations of De Marchi et al. (1996) sug- 
gest that this is the case, the above are hai^d hmits on p{0), due 
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to the assumption so far that Pint = 0. In fact Pint is likely to have 
some finite positive value, such that we have thus far underes- 
timated the accelerations felt by the pulsars. We now present 
some plausibility arguments that lead in general to more realis- 
tic estimates of a limit on p(0). 

Most millisecond pulsars known in the Galactic disk (1.5 < 
P < 9 ms) have inferred intrinsic surface dipole magnetic field 
strengths of 10*^- 10'' G, with an aver age of B 3 x 10** G (e.g., 
Camilo, Thorsett, & Kulkarni 1994j). For a given P and B, the 
period derivative is obtained (e.g., Manchester & Taylor 1977) 
from 

s -I ^ ■ ^ 

P\ / B 



P= 10" 



ms 



3.2 X lO^G 



(9) 



We have no a priori information about the magnetic field 
strengths of the pulsars in 47 Tuc, but if they were significantly 
higher than those for disk pulsars, so would the respective P. 
For a putative 6 ms pulsar, like 47 Tuc C, but with B = 10'" G, 
we obtain P = 2 x 10"'^ from equation ^ In turn, when used 
with equation |], this would imply p{Q) > 3 x 10^ pc"^, an 
unreasonably large value! A similar argument can be made for 
47 Tuc J. If the magnetic fields of the millisecond pulsars in 
47 Tuc are not likely to be much larger than those of millisecond 
pulsars in the disk, we may assume that they are, on average, 
comparable. 

Following the above arguments, we calculate from equa- 
tion|thelikelyvaluesofPintfor47TucCandJ: 1.7 x 10"^" and 
4.8 X 10"^° respectively, for B = 3.2 x 10^ G. These expected 
values can be decreased or increased by a factor of ten, consid- 
ering the entire range 10^ < B < 10'^ G. To be conservative, we 
consider only 10** < B < 3.2 x 10** G. For 47 Tuc C, the full con- 
tribution due to acceleration in the cluster (from eq. ||), is 5.1- 
6.7 X 10-2", while for 47 Tuc J it is 1. 4-5.7 x lO'^*'. These con- 
tributions are increased over the hard limits assuming Pint = 
by factors of 1.0-1.3 for 47 Tuc C and 1.5-6.1 for 47 Tuc J. 
Finally, the implied lower limits on p{Q) are increased to 3- 
4x10^MqPC-^^ and 2-lOx 10^ Mq pc-^ for 47 Tuc C and J 
respectively. 

We conclude from the present timing solutions that p(0) > 
3 X IO^MqPC"^. This lower limit can be compared to "stan- 
dard estimations" of p(0), ba sed on optical measurements of 
King model parameters (e.g., Webbink 1985). The aforemen- 
tioned revision in the core radius value for 47 Tuc implies 
that the existing estimate for /o(0), which scales as ( King 



1966; [Spitzer 1987 ), should be revised appropriately. Scal- 
ing the value for p(Q) tabulated by Webbink (1985), we find 
p(0) = 4 X lO^Mopc""*. While this is in agreement with our 
present lower limits, it is unlikely that both 47 Tuc C and J 
have the lowest magnetic field strengths of any known pulsars 
and so, in the context of a King model, the above discussion 
is strongly suggestive that for 47 Tuc, p{Q) > 4 x 10^ pc"-', 
and possibly is as much as a factor of two or so larger 

5.6. Pulsar Luminosities 

Except for 47 Tuc C, D, and J, pulsars in 47 Tuc are detected 
on fewer than 75% of observing days, and we have detected 
some of thepulsars on only one of 75 observing days (see Ta- 
bles ^ and ||). Given this fact, it is not appropriate to estimate 
average flux densities for most pulsars by simply averaging the 
detected signal-to-noise ratios and converting to flux densities 
using the approximately known system noise and detection ap- 
paratus characteristics. In this section we present a more rig- 
orous method to estimate mean flux densities of the pulsars in 



47 Tuc. 

The variability in apparent flux densities, apart from the oc- 
currenc e of eclipses, results from interstellar scintillation (R ick- 
ett 1977), and we proceed further by assuming that scintillation 
affects all pulsars in 47 Tuc similarly, since they are all located 
in relatively close proximity. 

Given this assumption, we construct cumulative histograms 
in \og{(f) for the frequently-detected pulsars in the cluster, and 
attempt to synthesize a "standard" histogram which we take to 
be representative of the relative flux variations due to interstel- 
lar scintillation for a pulsar in 47 Tuc. Initially the standard 
histogram is defined to be that of 47 Tuc J, the brightest pulsar 
The histograms of other frequently-appearing pulsars are then 
shifted by a scaling factor in a fitting procedure so as to match 
the standard as closely as possible, in a least-squares sense. A 
new standard is then defined where a of the «* strongest "ob- 
servation" is the average of the respective values from each of 
the individually shifted histograms, and this procedure is iter- 
ated until the scaling factors are stationary and a self-consistent 
solution is obtained. The scaling factors with respect to the fi- 
nal standard (see Fig. ^ therefore provide an estimate of the 
average relative signal-to-noise ratios of the various pulsars. 
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Fig. 6. — Cumulative "standard" histogram of \og{ait) for observations of 
pulsars in 47 Tuc. or is the signal-to-noise ratio of an "observation", (T,pScaled 
to the maximum signal-to-noise ratio, (Tmax, i.e., cfr = ajum'nx- See § p.6| for 
details. 



For 47 Tuc J we obtain the actual average a from the obser- 
vations (Table but correct it downwards to take into account 
the ^ 10% of observing days in which we do not detect the 
pulsar. We do this by assuming that on these days, on average, 
the pulsar has a flux density that is half of the lowest value de- 
tected, resulting in a 10% correction. After applying the scaling 
factors and converting to a flux density scale using the method- 
ology outlined in §^ we obtain the average flux densities listed 
in Table || for 14 pulsars. 

One obvious problem with this method is that for the less- 
frequently appearing pulsars only part of the cumulative his- 
togram is being fitted to — and to the extent that the standard 
histogram is in fact not precisely descriptive of the scintillation 
behavior of all pulsars, this may particularly skew the average a 
for such pulsai's. We have accounted for this in the uncertainties 
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listed in Table |[ In addition, until we calibrate the sensitivity of 
the observing system with a pulsed source of known strength, 
the conversion between cr and flux density is subject to a sys- 
tematic uncertainty that we judge to be of the order of ^ 25%. 

Pulsars 47 Tuc C, D, and J, for which we have the most reli- 
able flux density estimates, display maximum "daily" flux den- 
sities that are factors of 4-8 greater than the average values. 
This is the only information we have available that is of some 
use to estimate, very roughly, the average flux density for the 
pulsars that are detected very infrequently (Table ^). We obtain 
estimates of ~ 0.04 mJy for the average flux density of these 
pulsars. We note this plausible estimate simply to suggest that 
the flux densities of these pulsars need not be much lower than 
those of the weakest among the more frequently-appearing pul- 
sars. 

The weakest pulsars we have been able to detect in 47 Tuc 
have flux densities at 1400 MHz of about 5i4oo = 0.04 mJy, 
while the strongest, 47 Tuc J, has Sim) = 0.54 mJy. At the dis- 
tance of 47 Tuc (Table |l]), these correspond to luminosities at 
1400 MHz, Li4oo = Suood^, of 0.8-11 mJykpc^. Translated to 
luminosities at frequencies ^ 400 MHz, assuming a spectral in- 
dex for the pulsars in 47 Tuc of -2.0, these correspond roughly 
to 10 < L400 < BOmJykpc^. 

The limiting minimum luminosity for G alactic disk pulsars at 



400 MHz is abo ut 1 mJy kpc^ or lower ( |Lorimer et al. 1993 
Lyne et al. 1998[). For the Galactic disk ( Lyne, Manchester 



& Taylor 1985) and M15 ( [Anderson 1992^ the differential lu- 
minosity distribution for pulsars is found to be consistent with 
dN = L~' dlogL. Assuming that these findings apply to 47 Tuc, 
we estimate there are about 10 times as many potentially ob- 
servable pulsars in the luminosity decade 1-lOmJykpc^ as 
there are in 10-100 mJykpc^, where we already know of about 
20 pulsars. Hence we estimate a potentially observable popula- 



tion of pulsars in 47 Tuc, similar to the ones already detected, 
of about 200, above a minimum luminosity L400 ~ 1 mJykpc^. 
When beaming and selection effects are accounted for we judge 
it likely that 47 Tuc contains several hundred pulsars with 
^ ImJykpcl 

The actual determination of the luminosity distribution, and a 
detailed discussion of flux densities, spectral indices, and scin- 
tillation properties for the pulsars in 47 Tuc, will rely on more 
complete multi-wavelength data sets, and will be addressed 
elsewhere. It is already clear however, that we have detected 
only the few brightest pulsars in 47 Tuc, and continued searches 
of this cluster, particularly with greater sensitivity, should yield 
a bounty of new discoveries. 
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